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The uses of Whitham hierarchies 



Marcos Marinol*! 



Department of Physics, Yale University 
New Haven, CT 06520, USA. 



We review some of the uses of Whitham hierarchies in the context of the theory of the 
prepotential in M = 2 supersymmetric gauge theories. We focus on the structure of the 
contact terms in the twisted topological theory, and on the connection between Whitham 
hierarchies and the u-plane integrals for higher rank gauge groups, trying to put together 
the different approaches involved in this connection. We also review two other uses of 
the Whitham hierarchies: the interpretation of the slow times as supersymmetry breaking 
parameters, and the new techniques to extract instanton corrections using the RG equations 
written in terms of theta functions. 



§1. Introduction 



Shortly after Seiberg and Witten discovered the exact solution for the low energy 
effective action of M = 2 supersymmetric SU(2) gauge theories it was realized 
that this solution, as well as the generalizations ta other gauge groups, could be 
interpreted in the framework of integrable systemscP'EJ . This connection, which is 
| largely based on a common description in terms of Riemann surfaces and their Ja- 

cobians, has essentially involved two ingredients: the first ingredient has been the 
identification of an integrable classical mechanical system whose associated spectral 
curve reproduces the curve describing the low energy dynamics of the gauge theory. 



For example, for pure Yang-Mil]s theories, the relevant integrable system turns out 



to be the periodic Toda chaincf, while in the case n o: 
theory, the corresponding system is Calogero-Moser 



e mass deformed M = 4 
(which is in turn equiv- 



alent to the Hitchin system described in§)). The second ingredient is the theory 
of the prepotential: once the integrable system has been identified, one considers 
the quasiclassical Whitham hierarchy associated to the original hierarchy, which is 
constructed by introducing "slow" times instead of the original, "fast" times EP. The 
prepotential of the effective theory turns out to be, essentially, the logarithm of the 
quasiclassical tau function and hence depends on the slow times of the corresponding 



Whitham hierarchy 13 >E)>n 

Although there are general and rigorous arguments showing that the effective 
Af = 2 theories should be governed by integrable systemsLP'E9\ there is for the mo- 
ment no dynamical reason to explain why they are described by these particular one- 
dimensional mechanical systems. Nevertheless, this remarkable connection between 
two a priori unrelated fields has been very fruitful. For example, the connection to 
integrable systems gives a unifying approach to find the Seiberg- Witten curves for 
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the different gauge groups, because both Toda and Calogero-Moser systems can be 
formulated for any Lie algebra. In this paper, we will focus on the second aspect of 
the connection, namely the theory of the prepotential in the framework of Whitham 
hierarchies, and we will try to show that this framework is not only an elegant way 
to describe the prepotential, but also the appropriate point of view to understand 
some important aspects of M = 2 gauge theories. The main reason for this is that 
within the Whitham approach one can consider an "enlarged" prepotential with ex- 
tra parameters. Consider, for simplicity, the case of pure SU(N) super Yang-Mills 
theory. In this case, the usual Seiberg-Witten prepotential J 7 (a 1 , A) is a holomor- 
phic function of the coordinates of the moduli space, a 1 , where I = 1, ■ ■ ■ , N — 1, 
and the dynamically generated scale of the theory, A. The theory of Whitham hi- 
erarchies provides a generalized prepotential J 7 (a 1 , T n ), where T n , n = 1, • • • , N — 1, 
are the slow times of the Whitham hierarchy. This prepotential is a deformation of 
the Seiberg-Witten prepotential, in the sense that, if we put T\ = A, T n> \ = 0, we 
recover the original function J 7 (a 1 , A) . 

There are two contexts where deformations of the prepotential can be relevant. 
The first context is the study of contact terms in the twisted version of N = 2 gauge 
theories. These contact terms, as we will explain, appear when one computes the gen- 
erating functional of a certain class of operators, and the source terms for insertions 
of these operators can be regarded as deformations of the original action. The other 
context is the study of soft supersymmetry breaking terms, which can also be under- 
stood as deformations of the theory. In both cases, the deformations of the action 
can be described in terms of the Whitham time variables of the enlarged prepoten- 
tial. In this paper, we will mainly focus on the relation between the contact terms in 
twisted N = 2 theories and the theory of the prepotential from WJl i th a rn hierarc h ies 
The issue of contact terms has been previously addressed inO^llS'll^'llJ^^'ll^'llZP 
from different points of view. We will try to put together these approaches using the 
relation to integrable systems as a unifying framework. We will also review the uses 
of .Whitham hierarchies in the study of soft supersymmetry breaking, as explained 
in El. 

An extra bonus of the Whitham approach to the theory of the prepotential is a 
set of renormalization group (RG) equations for the prepotential in terms of theta 
functions associated to the root lattice of the gauge groupllJ. It was shown mx3 
that these equations give an elegant way to derive the instanton expansion of the 
prepotential following a recursive procedure. We will also review very briefly the 
strategy of the computation, and we will present some of the results obtained so far. 

The organization of the paper is as follows: in section 2, we explain the structure 
of contact terms in twisted N = 2 Yang-Mills theory. In section 3, we show how the 
contact terms can be understood from the point of view of Whitham hierarchies. In 
section 4, we study the connection between the blowup formula in twisted Ya ng- Mills 
theory and the Toda- Whitham approach. In section 5, we review, following Eif, the 
uses of Whitham hierachies in the problem of soft supersymmetry breaking in J\f = 2 
theories. In section 6 we present some results on instanton expansions using the RG 
equations. Finally, in section 7 we state some conclusions and open problems. 
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§2. Contact terms in twisted M=2 Yang-Mills theory 

2.1. Twisted J\f = 2 Yang-Mills theory 

One of the most important aspects of N = 2 Yang-Mills theories is their relation 
to Donaldson theory. It is a well-known fact that any N = 2 theory in four dimen- 
sions can be "twisted" to obtain a topological quantum field theory (TQFT), i.e. a 
quantum field theory whose correlation functions are formally metric- independent 
(for a review of TQFT and the twisting procedure, see for example E-P'Ef.) In the 
case of N = 2 Yang-Mills, the correlation functions of the twisted theory are in 
fact the famous Donaldson invariants of four-manifolds 15. In this section, we will 
consider some aspects of this twisted theory for the gauge group SU(N) on a simply- 
connected four- manifold X. In the following, r = N — 1 will denote the rank of the 
group. 

The first thing to do is to identify the gauge-invariant operators of the twisted 
theory. The twisted theory is characterized by the existence of a BRST charge 
or topological charge Q (which is in fact a particular combination of the original 
supersymmetric charges of the N = 2 theory), and the operators of the twisted 
theory have to be BRST invariant. These operators are called observables. The 
simplest observables of the theory are precisely the Casimirs of the gauge group, 
which we will take to be the elementary symmetric polynomials in the eigenvalues 
of the complex scalar field 4> of the N = 2 vector multiplet: 

Ok = -r Tv(j) k + lower order terms, k = 2, , , , ,N. (2-1) 

k 

Starting with these operators, one can generate the rest of the observables through 
the so-called descent procedure. To do this, one needs another operator in the twisted 
theory, the descent operator G^ which has spin one and comes from another com- 
bination of the SUSY charges. Acting with n times on the Casimir operators 
Ok one generates an n-form, The integration of this n-form on an ra-cycle on X is 
another observable, called the n-th topological descendant of Ok- As we are assum- 
ing that the manifold X is simply connected, the topological descendants will be 
constructed with two-cycles S in X: 

h(S) = I G 2 O k = \ I Tt(0*- 1 F) + . . . , (2.2) 
j s j s 

where F is the Yang-Mills field strength. The basic problem of the TQFT is to 
compute the generating function for the correlators of observables: 

Z(pk,fk,S) = (exp(j2(p k O k + fkh(S))p twisted theory • 

(2.3) 

k 

Notice that, in general, S will be an arbitrary linear combination of basic two-cycles 
Si, i = 1, . . . , b 2 (X), i. e. S = Yh=i ^ USi, therefore 



6 2 (X) 

h(S) = Uh(Si). 
i=i 



(2.4) 
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In total, we have r • 62 PO independent operators Ik(Si). 

The computation of ( |2.3| ) can be done using the low energy exact solution of 
N = 2. S 1 ?/ (iV) Yang-Mills theory. This solution is encoded in the hyperelliptic 
curve EHP 

y 2 = P 2 (x,u k )-4A 2N , (2.5) 

where P(x, u k ) = x N — Ylk=2 UkX N ~ k is the characteristic polynomial of SU(N), and 
Uk = {Ok) are the VEVs of the Casimir operators fl2.1|). This curve has genus g = r, 
and, as explained intf'cP, it can be identified with the spectral curve of the N site 
periodic Toda lattice. Associated to this curve, there is a meromorphic differential 
of the second kind, with a double pole at infinity, that can be explicitly written as: 

dSsw = P'(x,u k ) . (2.6) 

y 



This one-form satisfies the equation: 

ddSsw 
du k+1 



uJ k , (2.7) 



where to k = x k ~ l dx/y is a basis of holomorphic differentials. Let j 1 and 7d,j denote 
a symplectic basis of homology cycles for this curve, I = 1, . . . , r. The a 1 variables 
of the prepotential J- (a 1 , A), for the duality frame associated to the cycles 7 , are 
given by the integrals over these cycles of dSsw 



a\u k ,A) = ±-i - dx. (2i 



The same expression holds for the dual variables cldj with jd,i instead of 7 . Finally, 
the effective gauge couplings tjj are just the components of the Riemann matrix 
associated to the hyperelliptic curve (|2.5|). A fundamental aspect of the low-energy 
description of N = 2 gauge theories is that one has to choose a duality frame, which 
in the language of Riemann surfaces can be understood as a choice of the symplectic 
basis 7 , 7d j. Different duality frames are related by symplectic transformations of 
the form 

A B' 



F= {c D )tM2r,Z). (2.9) 

2.2. The integral over the Coulomb branch 

The moduli space of the hyperelliptic curve ( [2.5| ) is parametrized by the VEVs 
of the Casimirs u k . For some values of these VEVs, the curve will be singular. 
The singular locus is precisely the divisor D defined by the vanishing locus of the 
discriminant of fl2,5|) , A a = 0. It is well known that on this divisor there are extra 
BPS states that become massless. The Coulomb branch of the J\f = 2 theory is then 
defined as 

^Coulomb = C r - V. (2.10) 



The basic principle to compute the generating function ( |2.3|) has been introduced 
by Moore and Witten inEiP and states that this function is given by the sum of two 
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different contributions: one comes from the Coulomb branch, and the other comes 
from the divisor T>: 

Z = -^Coulomb + Z-D- (2-11) 

The first piece is given by an integral over the Coulomb branch, while the second 
piece involves, in general, Seiberg-Witten invariants associated to the extra massless 
BPS states. In this paper, we will focus on Zc ou iomb followingliHrEJ). As shown 



inEP'fiif, this contribution is different from zero only when b^X) = 1. The explicit 
expression of this integral can be worked out using the low-energy effective action 
encoded in (2.5) and reads 



^Coulomb =/ [dada]A(u k )XB(u k re^ u * +s2 Zf^ T *W. (2.12) 

J M. Coulomb 

The integrand of (E33) has various ingredients. First of all, tbe r e i s a gravitational 
part first studied inB' in the SU{2) case, and generalized in i3' t3> '113 to simply-laced 
groups. This part involves the factors: 

X/2 



A(u k ) 



aX ( det w) 

4 <r/8 



B{u k ) a = (3 a Ajl° . (2.13) 

The first factor is a modular form of weight (— x/2, 0), and B{u k ) is a modular 
invariant. In these equations, a, (5 are constants. Notice that the quantities involved 
here are very natural quantities associated to the hyperelliptic curve ( |2.5[ ), namely, 
the determinant of the matrix of periods of uj k , and the discriminant of the curve. 

Another ingredient in fl2.12|) is which is given by a sum over a lattice T. In 
order to write this quantity, let's introduce the following notation: 

du k 



z2fk da i- 

k 



(2.14) 



We also need some geometrical ingredients. As we are on a manifold of (X) = 1, 
given a metric g there exists a unique anti-self-dual form u £ H 2,+ (X,~R) such that 
to 2 = 1. The self-dual part of a two-form is then given by A+ = (A,w), where (•, •) is 
the usual product in cohomology. The lattice sum in & comes essentially from the 
evaluation of the partition function of the photons in the effective U (l) r gauge theory. 
A topological sector in the effective theory is specified by r 2-forms A G H 2 (X, R) . 
The lattice T consists of vectors of the form 



a = 



(2.15) 



i=i 



where {a/}/=i,..., r are the simple roots of S'[/(A^).0 In terms of these quantities, & 
can be written as follows: 

<f = (AetlmT)- 1,2 e^\—Vj\(linT)- 1 } JK V K Sl 
L 8tt 



*' For simplicity, we are assuming that there are no magnetic fluxes turned on. The general case 
is analyzed inE3. 



6 



Marcos Marino 



x ex P 
Xer 



-iirruiXi, X J + ) - iirruixL, X J _) - ivr(A • p, w 2 (X)) - iVj(S, XL) 



x / dr} 1 dx 1 exp< 
7=1 

1 



^^Wx J [47t(A? , u) + i(lmr) KL V L (S, u)] 



(2.16) 



where p is the Weyl vector, w 2 (X) is the second Stiefel- Whitney class of X, and rj , 
are Grassmannian coordinates which arise from the zero modes of the fermion 
fields in the theory. Tijk denote the third derivatives of the prepotential. Notice 
that \P depends explicitly on the metric of X through the two- form uj. It is precisely 
this dependence which leads to the wall-crossing phenomenaof Donaldson theory in 
manifolds with b^{X) = 1, when the gauge group is SU(2)t3. 

In ( |2.12| ) we have also included some terms of the form T^jS 2 in the exponential. 
These terms, which are proportional to the intersection form of the two-cycles, are the 
contact terms that we want to study. To understand the origin of these terms, notice 
that the operators Ik(S) that we have introduced in (|2.4D are non-local operators. 
These operators have a low-energy counterpart Ik(S) which can be obtained using 
the descent procedure in the effective abelian theory, and have been included in W '. 
But if we consider products of operators, and we go to the low-energy theory, we 
rather expect an extra contribution localised on the intersection of the surfaces that 
support the operators: 

ifcOSi)i fc OS 2 ) -» h(Si)i k (s 2 ) +r fc> ,(5i n s 2 ). (2.17) 

This is precisely what we have taken into account in (|2.12j) by introducing the contact 



terms T^/. It is important to notice that these terms are not predicted a priori by 
the Seiberg-Witten solution for the low-energy effective action. One has to find extra 
conditions in order to be able to find their structure. 

2.3. The structure of the contact terms 

To have more information on the contact terms, the first step is to write the 



integrand of ( 2.12j ) in a way that makes manifest the behavior under symplectic 
transformations. One introduces, then, the following generalized Siegel-Narain theta 
function: 



O r {Tu, a/, (3 1 ; P, £/) = exp 



-ivr(a 7 ,/3 7 ) + |^/,+ (Imr) /J ej, + - ^_(Imr) 7J ^-) 



x ex p 



-iTTTuCxLxi) - iTTTjj(X I _, X"L) - 2Tri(X I ,Z I ) + 2vri(A / ,a / ) 



(2.18) 



where X = X 1 + (3 1 . If we take 



2ir lbir 
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/3 7 = 0, ai = -W2(X), 1 = 1, •••,<• 
the lattice sum & can be written as 

C-2 

V = exp — Vi[CJmT)- x ] IJ Vj (det IuitY 1 / 2 

L87T J 



(2.19) 



IJ^Wexp -^ r ^ /Ji ,r ? / X J (Imr)^ L y L (5, W ) r (r 7J , a 7 , /3 7 ; P, fr). 

(2.20) 

It is easy to check© that W = exp — ^Vj[(lm.T)~ l ] IJ Vj \P is a modular form of 

weights (p2(X)/2, —1) (one has to take into account that b^(X) = 1). The Coulomb 
integral then reads: 



-^Coulomb = / [dudu] exp 

'-^Coulomb 



det(|^-)p(n fe )^(n fc )^. 



(2.21) 



The factor in the second line of (2.21) is a modular invariant (notice that, as X is 
simply connected, x = 2+62 (X))- We then see that this expression for the generating 
function (except for the exponential involving S 2 ) is the integral of a duality invariant 
object over a moduli space parametrized by the VEVs of the Casimirs, which are 
duality invariant coordinates. This is precisely what we expect on physical grounds: 
the generating function is a physical quantity and cannot depend on the choice of 
duality frame in the effective theory. This argument in fact forces the exponent in S 2 
to be a modular invariant as well, and this gives the first constraint on the contact 
terms: the quantity 



3fci + 



1 du k 



[(Imr) 



-1-1/J 



dui 



(2.22) 



87T da 1 da J 
must be invariant under the action of the symplectic group Sp(2r, Z). The other 
constraint on the contact terms has to do with their physical origin: these terms 
are quantum corrections and vanish at tree level, therefore they have to go to zero 
in the semiclassical region of moduli space (i.e. when A/a — ► 0). The duality 



transformation of the second term in (2.22) is easily worked out, and one finds that 
under a duality transformation it is shifted by: 



i duk 



{Ct + D)- 1 C 



1 J dm 



da J ' 



(2.23) 



The transformation of the contact term should compensate for this shift. Summa- 
rizing the discussion so far, we have found the following constraints for the contact 
terms: 

• They transform under an element of Sp(2r, Z) as follows: 



3fc / — * Tki + 



i duk 



(Ct + D)- l C 



1 J du 



da J 



(2.24) 
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• 2fc,z - ► semiclassically. 

As shown inEJ), these two properties determine the contact terms unambigu- 
ously. The problem, of course, is to find explicit expressions for them. 

In the case of the pure Yang- Mills theory with gauge group SU(2L the only 
contact term is T2 2- The solution to the constraints above was found intlP, and the 
result was expressed in terms of the Eisenstein series E^r). A general procedure 
to find Tip was described inllHP'E^P by using the so-called RG equ atio n of .hf = 2 
theories c3' c3 , as well as some aspects of the form alism developed int^'E^'Ef . The 
RG equation for J\f = 2 theories states thatEZP-Ef 

dT 1 , 

TfF„ = I" 2 ' (Z25 > 

where tq is defined as follows: for asymptotically free theories, A@ = e nlT ° , where (5 
is the coefficient of the one- loop beta function (for example, j3 = 2N for SU(N)). 
For the self-dual theories, tq is the microscopic coupling constant. The key point is 
to consider now the second derivatives of the prepotential. On the one hand, from 
( gjg ) one has that 

d 2 F _ 1 du 2 
aa J 0rt,~ Ada 1 ' [ > 

On the other hand, the fact that T^is invariant under duality transformations implies 
the following transformation lawt 



&T d^T nJK d^T 



*^ - d^ [(CT + D) 1 jC aw**- (2 - 27) 



This is precisely the shift we found in ( 2.24p under symplectic transformations. It is 



easy to see that d 2 J- /9tq is zero semiclassically (this is related to the fact that the 
classical prepotential is linear in the gauge coupling constant). It follows that the 
contact term associated to the quadratic Casimir can be written as 

T 2 ,2 = -^4- (2.28) 

T 2) 2 can be in fact evaluated very explicitly in the SU(2) case^'il). For the theories 
with Nf < 3, one finds: 

^=-^w(^) , +i(- +, "'4)- (2 - 29> 

and for the Nf = 4 case one has: 

T 2 ,2 = -^ E 2^)(^) 2 + Mro)l + \RE A (T ), (2.30) 

where R = J2a=i m a/2, and E 2 , £4 are the normalized Eisenstein series. This 
provides an elegant solution to our problem for this particular contact term. The 
authors ofEJ) found equivalent expressions for the SU(2) theories in the massless 
case and for Nf < 3. 
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The approach based on the RG equation suggests how to find Tki for any k,l. If 
one is able to find some additional variables T n , n = 2, . . . , N— 1 in the prepotential, 
which are invariant under duality and such that 

~n n+ i, (2.31) 



dT n 

then d 2 J- / dT^dTi will have, essentially, the duality properties of the contact term 
Tk+i,i+i, and one is halfway to the solution of the problem. The two obvious ques- 
tions are: Can we find an explicit construction of the generalized variables T n ?, 
and: Can we compute the second derivatives d 2 T '/ 'dT^dTi in terms of elementary 
data associated to the Seiberg-Witten solution? These two questions are answered 
inE-3) in the affirmative through the use of the Whitham hierarchy approach to the 
prepotential, which will be the subject of the next section. 

§3. Whitham hierarchies and contact terms 

3.1. Whitham hierarchy and Seiberg-Witten solution 

The approach to the^theory xtf the prepotential based on Whitham hierarchies 
has been developed inH^'B'B). IntJ it was shown that, using this approach, one can 



derive RG equations with the structure ( p. 31 ). We will follow heretB* and also®. 



The usual Seiberg-Witten differential is a meromorphic differential with a second- 
order pole at infinity, and such that its variations with respect to the moduli Uk are 
holomorphic differentials. In order to deform the Seiberg-Witten theory and to em- 
bed it in a larger framework, one considers a series of meromorphic differentials 
d£2 n (in the notation of I1ZP ) , with poles of order n + 1 at infinity and satisfying the 



condition 

ddft n 



holomorphic. (3-1) 



^moduli 

One then introduces a generating functional for these one-forms with auxiliary pa- 
rameters T n , n > 1: 

dS = ^T n df2 n . (3.2) 

n>l 

The parameters T n are precisely the slow times of the Whitham hierarchy. One of 
the results ofO) is an explicit expression for the meromorphic one-forms df2 n : 

A = /fW^, (3.3) 

where fllfcL-oo c k^j = c k^ k - It is easy to check that P+ N (x) = x, therefore 

dS(Ti,T n >2 = 0) = T\dSsw{A = 1). We then see that we recover the usual Seiberg- 
Witten differential when T n >2 = 0. Starting with the enlarged differential ( |3.2| ), one 
can construct a deformation of the usual Seiberg-Witten theory. For instance, one 
defines the periods of dS as 

a 1 (u k ,T x ,T 2 , ...) = ~^~j' I d ^ n 

ri>\ T 
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_ 7; r P{u k ) n fp'(u k ) 

T 1 a I (u k ,A = l) + 0(T n>1 ), 



dX 



(3.4) 



where a 7 (ufc,yl) is the usual Seiberg-Witten period of dSsw- In the same way, one 
defines the aoi with the expression (3.4) but with the jjj j replacing the 7 . Follow- 
ing now the usual steps in rigid special geometry, we obtain an enlarged prepotential 
^-(a^Tn), taking a 1 and the slow times as independent variables. 

We will now make more precise the relation with the usual Seiberg-Witten so- 
lution, followinglia) . We define the rescaled times T n and VEVs u k as folows: 



T n = T n /T?, u k = T?u k , 
with T\ = 1. It is easy to see that (3.4) can be written as 



a I (u k ,T 1 ,f 2 ,f 3 ,...) = V ^ / 

^ 2m J~i 

n>l ' 



Pl /N {u k )P'{u k ) 



P 2 (u k ) 



dX. 



4T 2N 



(3.5) 



(3.6) 



In particular, after setting T2 = T3 = 
a 1 {u k ,T x ,T n>1 = 0) 



... = we find that 
Tia 7 (ii fc ,l) = a 7 (n fc ,yl 



In conclusion, we may identify T\ with y1 in the submanifold T2 = T3 = ... = 0, 
provided the moduli space is parametrized with u k instead of u k . This shows that 
the Whitham hierarchy approach to the prepotential gives a deformation of the usual 
Seiberg-Witten theory, in such a way that the quantum scale can be identified with 
the first slow time of the Witham hierarchy. 

3.2. Derivatives of the prepotential and contact terms 

As we noticed in the last section, the slow times would be useful for the problem 
of the contact terms if they were "dual" to the higher Casimirs, in the sense of 
( |2.31| ). We should then compute the_derivatives of the prepotential J-(a,T). This 
computation was the main result oflHP. One of the basic ingredients in the answer 
is the Riemann theta function &[ct, P](£\t) and characteristics a = (a±, . . . ,a r ), 

P = (&,..., Pr)- 



0[aj](i\r)= exp[iTTT IJ (n I + (3 I )(nj + Pj) + 2TTi(n I + /3 I )^ I + ai / 
n K eZ 

The derivatives of the prepotential turn out to be: 

dT (3 



(3.7) 



dT n 2mn 



m+l,n+l ) 



d 2 T 

dT m dT n 



P dH n+1 
2mn da 1 

P 



(3.8) 



- — r 'Hm+l.n+l + 
2m 



p dn m+1 dn n+1 1 



ran 



da 1 da J 



ITT 



d TIJ \o g e E {o\T) 
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In these equations, Oe(0\t) denotes Riemann's theta function with a certain charac- 
teristic E, evaluated at the origin; 3 = 2N, m,n = 1, ...,r = N-r— 1, and derivatives 
with respect to T n are taken at constant a 1 . According toli^P , the characteristic 
E appearing in (3.10)-(3.11) is an even, half-integer characteristic associated to a 
particular splitting of the roots of the discriminant. Notice that the characteristic 
E depends on the duality frame. An explicit expression for E in the electric frame 
will be given in the next section. The functions l~im,n are certain homogeneous 
combinations of the Casimirs u^., given by 

Wm+i.n+1 = — resoo (p m l N {x)dPl IN (x)) = H n+1 , m+1 , 

and 

Um+i = T-im+1,2 = — reSooP™^ (x)dx = u m+ i + 0{u m ). 
m 

Here resp stands for the usual Cauchy residue at the point P. We have for instance 
^2,2 = 7~t2 = U2, ^3,2 = = U3 and 7^3 5 3 = U4 + ^if-u 2 ,. As we have seen, in 
order to recover the Seiberg-Witten solution, it is better to use the rescaled variables 
(A, T n >2). It is straightforward to reexpress all the formulae in (3.8) as derivatives 
of T with respect to a 1 ,T n and A. Most of the factors T\ can be used to promote 
Uk to iik or > rather, to the homogeneous combinations thereof: 

nj _ rpm+nnj . nj _ rpm+ln_l (q C\\ 

n. m +l,n+l — 1 1 itm+l,n+l => lim+l — J-\ n-m+l l^- y J 

(since 7i m +i = H m +l,2)- The remaining T±s are absorbed in making up a 1 = 
Tia 7 (ufc, 1) = a (uk,Tx). Altogether we find 

QT P sr t dT - I 3 sr t * 

dl^A ~ 2ri A* mT mT n n m+1>n+1 , mT m H m+1 , n+1 , 

° n,m>l ^^n m>l 

d 2 T _J3_s rf dttm+i d 2 T _ 3 dH n+1 

^ = -77- 2^ TmTn—TTZJ TTTj- — 9 Tl j log 9 E (0 \t) , (3.10) 

a loErylr 27tz aa J oa J vk 

y ° ' m,n>l 

d 2 T b 2 dn n+1 _ . a# m+1 i 

L Tm ^j- log <S>^(0|^) , 



<91og/ldT n 2ttot «9a J ^ 5a J 

d 2 T (3 (~ B dH m+1 dHn+i 1 A 

with m,n > 2. In these expressions T\ = 1. The restriction to the submanifold 
T2 = T3 = ... = yields formulae which are expressed in terms of the original 
Seiberg-Witten data. Notice that in this subspace a 7 (ufc, T\, T n> \ = 0) = a 1 , hence 

9T 3 ^ dT 3 ~ 

rn-2 , —5- — - — —rt n+ i , 



91ogyl 2iri ' df n 2mn 
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d 2 T 13 dii 2 d 2 T _ (3 dH n+1 

dtfdlogA ~ 2m da 1 ' d&dfn ~ 2mn da 1 

d 2 T p 2 dH 2 dH 2 1 



d(log A) 2 2m da 1 da J m 
d 2 T f3 2 dH 2 dH n+1 1 



d TlJ \ogO E (0\r) , (3.11) 
d Ti Aog0 E (O\r) , 



dlogAdf n 2mn da 1 da J in 

d 2 F (3 P dH m+1 &H n+1 1 \ 

~ = -7T - Wm+l,n+l H Jw — Orj j log6>B(0 T) 

8T m dT n 2m \ mn da 1 da J m J 

with m,n > 2. We will always work with the coordinates 7l and T n defined in 
and, therefore, hats will be omitted everywhere. 

The first derivatives of the prepotential with respect to the slow times can be 
interpreted as generalized RG equations. These equations show that, indeed, the 
slow times T n are dual to the Casimirs Ti n +l, at least in the submanifold T n > 2 = 0. 
One can also work out the transformation properties under the symplectic group 
of the second derivatives of the prepotential. To do this, notice that the quantities 
TCfi^m are combinations of gauge invariant VEVs, therefore they are invariant under 
duality. The slow times are deformation parameters of the theory, and they are 
also invariant. The symplectic transformation of the periods a 1 , ao,i are the usual 
ones for the variables a 1 , cld,i of the Seiberg-Witten ^eory, and one derives the 
transformation law for the prepotential obtained int 



.27 



T r = T + ^a I (D T B) ij a J + \a DI (C T A) IJ a D ,j + a I (B T C) I J a D>J . (3.12) 



One then finds 



0>. 



■IVr + DW" (3,3) 



8T n dT m 8T n dT m 8T n da lVK ' iJ dT m da K ' 

This transformation law can also be checked on the explicit expressions in (3.10), by 
using the transformation law of the theta function under a symplectic transformation: 

0[a, 0} {£\t) -» e^(det(Cr + D)) 1 ' 2 exp Un? (Ct + D) _1 C^l@[a, P]{£\t), (3.14) 

where is an eighth root of unity. The transformation law ( |3.13j ) is the gener- 
alization of ( |2.27| ) that we were looking for. We would be tempted to identify the 
contact term Tfc+i f + i with the second derivative d 2 J-/dTf t Ti, up to a normalization 
factor. There are, however, two subtleties one has to take into account in order to 
make this identification. The first one is that, in the above model, the Casimirs are 
not Uk, but the homogeneous combinations TC^- Therefore, to identify the contact 
terms with the second derivatives of the prepotential, we have to define the operators 
Ik(S) starting with TLk instead of the Uk in the descent procedure. This is simply 
a redefinition of which gauge-invariant operators we take as a basis to construct 
the observables. The second subtlety is that the contact terms are not only deter- 
mined by their transformation properties. As we discussed in the previous section, 



The uses of Whitham hierarchies 



13 



there is the extra constraint that T^i vanishes semiclassically. Let's then analyze the 
semiclassical behavior of d 2 J- '/dT^Ti. If we use the semiclassical expression for the 
effective gauge coupling 77 j, we find that: 



d TlJ lo g e E (0\T)~O[ , (3.15) 



where Z$ + = a - dq_, o = J2i aI ®ii an d c?+ are the positive roots. Therefore, the 
term in d 2 J- '/dT^Ti involving this derivative will vanish in the semiclassical region 
(A/Z$ + ) —* 0, for k, I = 1, . . . , N — 1 (we take T\ = A here). On the other hand, 
Wn+i,m+i does not vanish. Therefore, d 2 T ' /dT^Ti does not have the required behav- 
ior for k, I = 2, . . . , N — 1. This can be cured by a redefinition of the prepotential. 
Consider the "reduced" prepotential 

■pied _ jr _|_ ^2 T n T m H n +l,m+l- (3.16) 

n,m>2 

The second derivatives of this prepotential with respect to the slow times are, in the 
Seiberg-Witten submanifold T n >2 = 0, 

9 2 ^ red \ /T dH n+ i 

d TlJ logO E {0\T), (3.17) 



<dT n dT m /T n > 2 =o 2mn da 1 da J m TlJ 
and the contact terms are given by 

_ Trikl ( d 2 F cd \ 
T k+1 , l+1 - _^^j T ^ a=o . (3.18) 

Notice that the piece that we substracted in the prepotential has a nonzero classical 



limit, while (3.17) is a quantum correction. In other words, the classical limit of the 
original prepotential T includes deformations which are quadratic in the slow times. 
This is not what we expect in the problem of evaluating the generating function 
( |2.3| ). In this generating function, all the deformations are linear in the sources, and 
the quadratic terms in these sources that appear in the effective theory (the contact 
terms) are only due to quantum effects. 

To summarize, we have seen that the prepotential defined in terms of the 
Whitham hierarchy is the appropriate object to understand the contact terms of 
the twisted theory. However, the above remarks suggest that there should be a way 
to improve the theory of the prepotential, in such a way that one obtains the reduced 
version introduced in ( p. 16 ), and that the slow times T n can be associated to any 



homogeneous combination of the Casimirs of degree n + 1. 

§4. Blowup formula 

In this section, we will discuss a further connection between integrable systems 
and the Coulomb integral, which appears when one considers the blowup formula. 
There are three reasons why such a formula is interesting: first of all, as shown inl!3 , 
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the_behavior under blowup can_hfi used to derive the contact terms. Second, as shown 
mi3 and further discussed inll^, the blowup formula involves in a direct way the 
tau function of the periodic Toda lattice which underlies the Seiberg-Witten curve. 
Finally, using the blowup formula one can fix the characteristic E which appears in 
(3.8) Ej). Of course, this formula is also very important in the context of Donaldson 
theory, as it makes possible to relate the Donaldson invariants of a manifold and its 
blowup 0. This formula was obtained using the u-plane integral in0\ for the case 
of SU(2), and generalized to the higher rank case inlij). 

Suppose that we have a four-manifold X, and we consider the blownup manifold 
at a point p, X = B1 P (X). Under this operation, the homology changes as follows!^: 

H 2 (X) -» H 2 (X) = H 2 (X) ®Z-B, (4.1) 

where B, the class of the exceptional divisor, satisfies B 2 = —1. As b 2 (X) = b 2 (X), 
it makes sense to ask how the intergrand of ^coulomb changes under the blowup. 
First of all, the blownup manifold X has an extra two-homology class, therefore 
there are extra operators Ik(B) that must be included in the generating function. 
We will then write S = S + tB. We want to compute then: 

Lxp[j2{fkIk(S)+t k I k (B)+p k G k )]\ , (4.2) 

\ k IX 

in terms of correlation functions of the twisted theory on X. In this equation, 
tk = fkt- The first thing to do is to analyze the change in the integrand under blowup. 
This is easy to obtain if we also assume that the metric is such that (B,u>) = 0. In 
this case, the lattice sum changes as follows: 



*x = I ^exp^ir/jnV + i^^-^jn 1 -ivr^Vj W x , (4.3) 

V n 1 k I / 

and we see that 

9 jt = 9[£,0](&)Vx, (4.4) 
where 0[A, 0](£|r) is a theta function with 

& = E|£t£. (1/2,..., 1/2). (4.5) 

k 

Notice that we have extracted from the Siegel-Narain theta function a standard 
theta function on the hyperelliptic curve ( |2.5| ), The above expression is valid in the 
electric frame, and the characteristic of the theta function is inherited from the term 
(p-X,w 2 (X)) in (2.16). 

Let's now analyze the measure in the integrand. As x(A) = x(A) + 1 and 
o~(X) = o~(X) — 1, the measure picks an extra factor 

(detpL) 1/2 A -Vs = ^l , (4 .6) 



*-* It is interesting to notice that the derivation of this formula inH- 1 was one of the first hints 
of a relation between Donaldson theory and elliptic curves. 
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as a consequence of Thomae formulae. Putting all these factors together, we see that 
the blowup factor in the integrand is given (up to a constant) by: 

rfey = e-EWM gM (47) 
[k]k) 0[AO](O|r)' ( j 

As (|4.7|) is an extra factor in the integrand of ^coulomb) it follows by our arguments 
above that it must be invariant under duality transformations. Using (3.14), one 



can easily prove that the duality invariance of the blowup factor fixes the duality 



transformation of the contact terms, and one finds ( 2.24 ) again. The reason for the 
notation in ([4.7]) is that the blowup factor is essentially the tau function of the Toda 
hierarchy (see, for example, Notice that, in this identification, the coupling 
constants t k are interpreted as the fast times of the hierarchy. Recall that there are 
two points of contact between integrable systems and the Seiberg-Witten solution: 
one is the fact that the hyperelliptic curves are spectral curves of periodic Toda 
chains (for pure Yang-Mills), and the other is that the prepotential is the logarithm 
of the quasiclassical tau function. Now, we see that the tau function itself appears 
in the twisted version of the theory as a blowup factor. This also gives a very direct 
link between Donaldson theory and its generalizations, and the hyperelliptic curve 
encoding the Seiberg-Witten solution. 

In principle, the blowup factor depends on the moduli of the hyperelliptic curve 
in a complicated way. However, it turns out that the blowup factor only depends 
on the fast times % and on the Casimirs, as we have indicated with our notation. 
More precisely, let ft = (ri2, • • • ,un) be a vector of nonnegative integers, with |n| = 
ri2 + • • • +un- Then, there are polynomials Bfi(u2, ■ ■ ■ , un) in the Casimirs such that: 

r(t k \u k )= J2 t 2 2 ---tN N M^2,...,u N ). (4.8) 

|n|>0 

The physical reason for this, as explained in0>,is that one can interpret the blowup 
as a local "defect" created by an analog of the puncture operator. The effect of this 
puncture associated to the exceptional divisor should be represented by an infinite 
number of local observables. The ring of local observables is in fact generated by 
U2, ■ ■ ■ ,un, and an identity like (|4,8| ) should be valid. One can_Lpdeed prove (|4.8j) 
by using duality invariance, i?-charge arguments and regularity^. A consequence 
of (Wjfy is the following relation between the generating functions: 



Uxp[£(fkIk(S)+t k I k (B)+ Pk O k )] ) = Uxp[£(fkIk(S)+PkO k )\T(t k \O k )) . 

\ k I X \ k IX 

(4-9) 

Notice that the blowup factor r(t k \u k ) does not depend on the manifold X. We have 
then proved that the blowup changes the generating function by a universal factor 
which can be expressed in terms of the zero observables. 

As a corollary of (fi]^), we can derive an explicit expression for the contact terms. 
The fact that the behavior under blowup can be used to obtain the contact terms 
was first remarked in©. To obtain the expression, we simply expand (4.7) to second 
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order in tk- The first derivative of the theta function is zero due to the choice of 
characteristic. We then find: 

r(t k \u k ) = 1 - Wt m + ±-.d TlJ l og 0[4 01(0^)^0]^ + • • • , (4.10) 

k,l \ / 

Because of (4.8), this means that 

T k ,i = -^-dru log 0[A, O](O|r)^0 + S^> 2 , . . . , u N ), (4.11) 

where n^i are the vectors with \fik,i\ = 2, corresponding to the quadratic terms in 
(|4.8| ). The requirement that Tj^/ vanishes semiclassically implies that B^ k ,(w2, • • • , un) 
0, and we finally find: 

r t ,, = -^,,lo ge[ AOl(O|r)^0, (4.12) 

which is the expression found in©. If we compare now to ( |3.18| ), we see that they 
have the same structure. The only difference is due to the fact that the expression in 
( 3.18| ) is valid if the descent operators are constructed with the operators TLk- Here, 
instead, we are considering the operators that we defined in (|2.2|) starting from O k - 
An interesting consequence of this rederivation is that we can read the characteristic 
E in (3.8) from the expression (4.12): 

E = (A,0). (4.13) 

We will provide another check of this identification in the section 6. 

To clarify our statement about the existence of the polynomials in Q4,8| ), it is 
useful to consider in some detail the SU(2) case. The blowup factor is simply given 
by0>: 

rm=e -t^^^l } (414) 

where $4 is the Jacobi theta form with characteristic [1/2,0], and we put u = u<i- 
The quotient of these theta functions_£an be written in terms of a Weierstrass sigma 
function using the following identityc2P : 



a 3 (oj 2 z), (4.15) 



04(O|t) 

where 0J2 is the a-period (i.e., r = U1/0J2), and 772 = ((U2/2). Using the expression 
for the contact term in ( 2.29| ), the identity 

m = t—E 2 (t), (4.16) 

and the fact that u>2 = (8-ir/\f2)(da/du), one finds: 

r(t\u) = e-* 2 tcT3(i|) (4.17) 
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This is precisely the expression found in© in the context of Donaldson theory (for 
zero magnetic flux). The key fact is that the sigma functions &i(t) can be expanded 
around the origin, and the coefficients of the Taylor expansion are polynomials in the 
root ej, and the functions g2 , 53 (in the Weierstrass description). These quantities 
depend only on u. From the Seiberg-Witten curve for the SU(2) pure Yang-Mills 
theory one finds in fact: 

U (4.18) 



63 12' 92 4V 3 4/' 93 48 V 9 4> 
Therefore, the expansion of ( |4.17D is indeed of the form ([4.8]), and using that 

a 3 (t) = l-e 3 t 2 + 0(t*), (4.19) 



it follows that the quadratic term in t in ( 4.17] ) is zero, as expected from ( |4.11| ) 



If we consider now the higher rank case, the above remarks on SU{2) suggest 



that the expansion in (4.8) should involve some kind of hyperelliptic generalization 
of the sigma functions. It is also quite possible that the interpretation in terms of the 
Toda- Whitham hierarchy gives a constructive way of computing these polynomials. 



§5. Soft Super symmetry Breaking with Higher Casimir Operators 

In this section, we will give a very rough overview of another use of Whitham 
hierarchies: soft supersymmetry breaking witlLspurion superfields. 

The spurion formalism was introduced inH and it is a very useful tool to break 
supersymmetry in an explicit way. The starting point of the spurion formalism 
is a coupling constant in a supersymmetric Lagrangian, call it m. This coupling 
constant can be promoted to a superfield, m — * M, and this will give another 
Lagrangian which will be supersymmetric as well. The superfield M is called a 
spurion superfield. Notice that the original Lagrangian is recovered by taking the 
VEV of the scalar component of M to be equal to m, and setting the rest of the 
fields to zero. To break supersymmetry, one gives a VEV to an auxiliary field in 
M. The resulting Lagrangian will be nonsupersymmetric due to the extra terms 
generated in this way. An interesting example is the mass term for the N = 2 quark 
hypermultiplet, which in N = 1 superspace has the form: 

J d 2 9mQQ, (5.1) 

where Q, Q are the two Af = 1 chiral superfields that correspond to the N = 2 quark 
hypermultiplet. If we promote m to an M = 1 chiral superfield: 

m ^ M = (j) m + V29ip m + 9 2 F m , (5.2) 

we obtain an J\f = 1 Lagrangian, and we recover the original one if we put {(j> m ) = m. 
If we want to break supersymmetry, we take: 



M = m + 6 2 F 1 



(5.3) 
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This induces an extra term F m qq (a mass term for the squark) which breaks super- 
symmetry down to ftf = 0. The VEV of the auxiliary field, F m , becomes a SUSY 
breaking parameter. Notice that one can also take m = 0, F m ^ to generate a mass 
term for the squarks while keeping the quarks massless. This is the usual procedure 
to decouple the squarks through soft supersymmetry breaking. 

We can try to proceed in the same way with M = 2 gauge theories described by 
a prepotential. In the most general case, the prepotential is a function 

F = F(a I ,m f ,T n ), (5.4) 

where T n are the Whitham slow times, and T\ = A. In this Lagrangian, the variables 
a 1 correspond to N = 2 vector superfields, while the rest of the variables are cou- 
plings. Any of these couplings can be promoted to an J\f = 2 vector superfield, and 
then one breaks super symmetry b y g iving VEVs to the auxiliary components. This 
strategy was followed in^)>E«M> in the context of the original Seiberg-Witten 
theory, when T n = for n > 2. The only coupling constants are then A and rut. In 
the enlarged prepotential, we can construct a more general, nonsupersymmetric de- 
formation of Seiberg-Witten theory, by promoting T n >2 to spurion superfields while 
setting the VEVs of their scalar components to zero (for T\, the scalar component 
has to be different from zero and equal to the quantum scale A). It is convenient to 
define the couplings s n as 

sx = -ilogA, s n = -iT n , n = 2,...,r. (5.5) 

We then promote these couplings to J\f = 2 vector superfields S n . Such a superfield 
can be written, in J\f = 1 language, as one M = 1 chiral superfield (that we will also 
denote by S n ) and one M = 1 vector superfield V n . We then have, in terms of J\f = 1 
superfields, 

Si = si + fiFi , V, = Vi = ^D^-P , (5.6) 

8. = <PF. , V n = \DJH> , » > 2. (5.7) 

The Whitham times are then interpreted as supersymmetry breaking parameters, 
and supersymmetry is broken down to Af = 0. To obtain the explicit form of the 
soft breaking terms, we can expand the prepotential around S n >2 = 0. The terms 
of order 0(S 3 ) in this expansion do not give any contribution to the Lagrangian, 
because they involve too many 6's and the integral in superspace will vanish. To 
analyze the softly broken theory, it is enough to consider the terms which are at 
most quadratic in the slow times, and the following expression is exact: 

n>2 - n,m>2 - 

(5-8) 

The first term in this expansion is just the Seiberg-Witten prepotential, after setting 
T\ = A. Notice that the prepotential defined in the Whitham framework is defined for 
the effective theory. Nevertheless, it is easy to see what is the classical prepotential by 
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going to the semiclassical region and switching off the quantum corrections. Again, 
it is more convenient to use the reduced prepotential ( 3.16| ). In this case, the second 
derivatives appearing in (|5.8|) vanish semiclassically, as we showed in ( |3.17l ). The first 
derivatives can be read from (3.11). They are written in terms of Casimir operators 
and therefore have a well-defined classical limit. We finally obtain the following 
expression for the reduced, classical prepotential: 



1 



•r-red _ 
class 27T ^ n 

n=l 



S n TL 



n+1, 



(5.9) 



where the spurions S n are given in fl5.6Q — (5.7). The microscopic Lagrangian asso- 
ciated to ( ]5.9| ) can be written as follows. First, one defines a generalized matrix of 
couplings as follows: 



(class) 

T ab 



f)2 •Tied 
u ^class 
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q2 •rrcd 
u ^class 



(class) 



o , 



(5.11) 



(5.12) 



dS m dS n 

where the dots in eq.(5.11) denote the derivative with respect to (f) a of lower or- 
der Casimir operators. The indices a, b, c, ... belong to the adjoint representation of 
SU(N), and are raised and lowered with the invariant metric, and n, m, ... corre- 
spond to the variables s n . The complex scalar field of the N = 2 vector multiplet is 
written as <p = J2 a 4> a T a , where T a is a basis for the Lie algebra of SU (N). Finally, r 
is the classical gauge coupling, and it is related to the spurion s\ = —i log A through 
the one loop formula A@ = e nlT . We define now: 



lass) 



4tt 



-Im r 



(class) 



We find 



I Q T (dass) r. 

— Im — 

8vr d6 a 



(r^ b )F m + (\ a X b )F m + iV2(X a ^ b )D r , 



where BJ"^^ is the classical value of the duality invariant quantity 



(5.13) 



C = C N=2 - B^l) UmF: + \D m D n \ + f% c b^ a V^D m <j> b p 



(5.14) 



B mn = b m b ~lab b n _ h mn 



(5.15) 



and f a bc are the structure constants of the Lie algebra. In (5.14), A, V are the 
gluinos and <f> is the scalar component of the N = 2 vector superfield. We see that 
the spurion that corresponds to s\ gives mass to the gauginos of the J\f = 2 vector 
multiplet and to the imaginary part of the Higgs field eft. The spurions corresponding 
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to higher Casimirs, on the other hand, give couplings between the Higgs field and 
the gauginos. Notice that the spurion superfields S n have dimensions 1— n, therefore 
the supersymmetry breaking parameters F n , D n have dimension 2 — n. For n > 2, 
they will give nonrenormalizable {i.e. irrelevant) interactions in the microscopic 
Lagrangian. This does not mean that the resulting perturbations do not change 
the low-energy structure of the theory: the operators we are considering can be 
dangerously irrelevant operators, as in the related theory analyzed inEj), and in this 
case they will affect the infrared physics. 

One of the advantages of this procedure to break supersymmetry is precisely 
that one can also write the exact low energy effective theory associated to (5.14), up 
to two derivatives. This is due to the fact that the dependence of the Lagrangian on 
the SUSY breaking parameters is controlled by the dependence of the prepotential 
on the spurions, and this is also given by the Seiberg-Witten solution. Notice that, in 
the low energy theory, the couplings ( |5.10 )-( p.l2 ) receive quantum corrections, and 
can be computed in terms of the hyperelliptic curve data. Using these couplings, 
one can write down an effective potential and study the vacuum structure of the 
theory. The effect of the nonsupersymmetric deformations associated to the higher 
order Casimirs in the vacuum structure of the theory has been explored int^P. 



§6. Instanton Corrections 

In this section we will briefly review another use of the Whitham hierarchies: 
the computation of instanton corrections. For more details, seec3\t3^ andfcf. 

One of the main results of the Whitham approach to the theory of the prepoten- 
tial have been the equations for the first and second derivatives of the prepotential 
in (3.8) and (3.10) derived in0. InEzI it was realized that the RG equation 

dlogyl 4m 

is very useful to derive a recursion relation for the instanton contributions. In or- 



der to compute the instanton corrections, however, (6T) is not sufficient and one 
needs additional input. This is usually provided by the Picard-Fuchs equations for 
the periods. The Picard-Fuchs equations are difficult to derive and solve when the 
rank of the gauge group is larger than one, although techniques from topological 
Landau-Ginzburg theories can make them more instrumental in order to obtain the 
one-instanton correction to the prepotential for the ADE series©. It turns out that 
the equation for d 2 F/d(log A) 2 in (3.11), together with (|6.1|) , provides enough in- 
formation to obtain the instanton expansion of the prepotential in the semiclassical 
region to any order, and we don't have to make use of the Picard-Fuchs equations. 
As we will see, the connection of SU(N), N = 2 super Yang-Mills theory with 
Toda- Whitham hierarchies embodies in a natural way a recursive procedure to com- 
pute all instanton corrections. The essential ingredient that makes this possible is 



The uses of Whitham hierarchies 



21 



the relation of the derivatives of the prepotential with the theta function associated 
to the root lattice of the gauge group. To see how this works, consider the instanton 
expansion of the prepotential: 

* = ^0 E Zl + J. £ Zi+ log ^ + ^ £ (0.2) 

/c=l 

In this expression, 2~)a + denotes a sum over positive roots. The expansion is in 
powers of A@, where f3 = 2N for SU(N), and k is the instanton number. We then 
have 

^ = 1 Lg (Wi (Z M - (6.3) 

which, according to (3.11), should be equal to 
d 2 F d 2 dU 2 dU 2 1 

The derivative of the quadratic Casimir has also an expansion that can be obtained 
from the RG equation and ( |6.2| ): 

C/ja J + 53*^,7^* (6-5) 



da' /? drfdlogA 



k=l 



where Cjj is the Cartan matrix, J-^j = dJ-^/da 1 , and we have taken into account 
that ^ J2s+ = \o J I Cija J . The couplings in the semiclassical region are ob- 
tained again from the expansion ( |6.2[) : 



^ = ^E^r^riogl-^l+ — ■ (6.6) 

a + \ / k=l 

with J^/j = d f^j . For convenience, in a term ^ -^r~^T (iT 7 " ~ 3 
has been set to zero by a suitable adjustment of the bare coupling 2-7riro = 3N. If 
we set a = 1 f^ai and define 

a-T'l-a = n^kjjn 3 , (6.7) 
1, J 

we see that the theta function Oe in the semiclassical region can be written as 

00 00 / 00 -. \ 

&e(o\t) =ee (-!)** n<* 5+) n e si ^'-^ ^ fcm u 27vr - 

r=0ae4 r (5+ fc=l \m=0 ' / 

(6-8) 

In the previous expression, /> is again the Weyl vector, and A r C A is a subset of the 
root consisting of the lattice vectors a that fulfill the constraint J2d + (5-a+) 2 = 2Nr. 
In particular Z\i is the root system, i.e. the simple roots together with their Weyl 
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reflections. In the above expression, we have used the characteristic ( |4. 13 ) , which 
is the appropriate one if we are working in the electric frame, as we should do in 
the semiclassical region. If we insert the expressions ( |6.5[) and ( |6.8[) into (^J), and 
we equate the coefficients of the different powers of A, we find a set of recursive 
equations for T k which make possible to compute all the instanton coefficients by 
starting from the perturbative contribution to the prepotential. For example, one 
finds for the one-instanton correction: 



(a-a+) 2 



(6.9) 



and for the two-instanton contribution: 



?1 



E(-i)^n^T 



T x + 2(a-F' k )Z d + -{a-F'{-a)Zl 



E 

/3ez\ 2 



(P-d+f 



(6.10) 



where a(F' k 



J2i nI J~kj- The above expression makes patent the recursive character 
of the procedure. We then see that, with this method, one can find the instan- 
ton corrections to the prepotential in a very simple way The explicit expressions 
Q6.9| )-(6.10) are in full agreement with the results offiZP't^'ElP (seet3^ for a detailed 
presentation). This agreement gives a further check of the RG equation (3.11) and 
of the choice of the characteristic ( 4.13] ) . 

These results on the instanton corrections can be extended in many ways. The 
equation (3.11) can be also used to study the strong coupling regime near the points 
of maximal singularity^. Moreover, one can analyze with this technique all the 
classical groups with or without matter content to find general expressions for the 
instanton corrections US. 



§7. Conclusions and Outlook 



We have seen that the approach to the N = 2 prepotential based on the theory 
of Whitham hierarchies gives a very useful framework to understand deformations of 
N = 2 supersymmetric gauge theories. These deformations arise in different physical 
contexts. We have explored in detail the structure of the contact terms in the twisted 
version of J\f = 2 theories, and also the nonsupersymmetric deformations associated 
to soft supersymmetry breaking. In both cases, the approach based on Whitham 
hierarchies provides the right conceptual framework and the technical tools to derive 
the precise form of the deformed theory. We have also seen that one can obtain the 
instanton expansion in the semiclassical region by using the equations for the first 
and second derivatives of the prepotential with respect to the quantum scale. 

Nevertheless, there are many issues that should be further clarified. The frame- 
work of Whitham hierarchies should be generalized to other g a u g e groups and matter 
content. Some steps in this direction can be found in In the context of 
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the twisted theory, one would like to introduce slow times associated to any homo- 
geneous combination of the Casimirs, and construct a prepotential depending on 
these slow times and with the right semiclassical behavior. In this sense, although 
the explicit construction given mx3 has clarified the theory of the prepotential and 
its applications to Donaldson theory, one should be able to improve it along the 
lines that we have suggested. It would be also nice to have an a priori connection 
between the Whitham hierarchy approach to the prepotential, and the structure of 
the generating function for the twisted theory. We have followed a rather indirect 
approach, based on a set of constraints for the contact terms and the RG equations. 
The approach of©, based on Hamiltonian deformations of the prepotential, should 
be useful in deriving this connection. 

It would be also very interesting to find an explicit construction of the polyno- 
mials appearing in ([|]^) for the higher rank case. As we have suggested, this con- 
struction will involve interesting generalizations of sigma functions for hyperelliptic 
curves. Another reason for addressing this issue is that the structure of the blowup 
formula gives a very direct connection between the Toda hierarchy underlying the 
Seiberg-Witten solution, and the generalizations of Donaldson theory. This relation 
is certainly intriguing and suggests that many structures that have been found in 
two-dimensional topological gravity and two-dimensional topological matter could 
be also relevant in four dimensions. 
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